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We estimate the long-distance contribution to the width difference of Bs—Bs system, based mainly 

on two-body Da Di* modes and three-body Di* D^*' K^*' modes (and their CP conjugates). Some 

higher cs resonances are also considered. The contribution to AFs/Fs by two-body modes is (10. 2± 

3.0)%, slightly smaller than the short-distance result of (13.3 ± 3.2)%. The contribution to AFs/Fs 

by D*o(2317), Dsi(2460), and Dsi(2536) resonances is negligible. For the three-body oi*^ D'-*^ K'-*^ 

^^ . modes, we adopt the factorization formalism and model the form factors with off-shell Di* poles, 

Cn ' the _Dsj(2700) resonance, and non-resonant (NR) contributions. These three-body modes can arise 

{3J[), through current-produced or transition diagrams, but only SU(3)-related D*]j^D^*'K modes from 

^ ' current diagram have been measured so far. The pole model results for D^^^D^*^ K agree well 

with data, while Du,dD^ K rates agree with data only within a factor of 2 to 3. All measured 

Qs ' D^^D^*'' K rates can be reproduced by including NR contribution. The total AF^/Fs obtained is 

(16.7 ± 8.5)%, which agrees with the short-distance result within uncertainties. For illustration, 

we also demonstrate the effect of _Dsj (2700) in modes with D^*'K*. In all scenarios, the total 

AFs/Fs remain consistence to the short-distance result. Our result indicates that (a) the operator 

product expansion (OPE) in short-distance picture is a valid assumption, (b) approximating the 

^L(' Bs — >■ Ds* Da* decays to saturate AFs has a large correction, (c) the effect of three-body modes 

QJ , cannot be neglected, and (d) in addition to Da and D*a poles, the Dsj (2700) resonance also plays an 

i-C^i important role in three-body modes. Future experiments are necessary to improve the estimation 

of AFs from long-distance picture. 
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I. INTRODUCTION AND MOTIVATION 

One of the most exciting news in_particle physics last year is the anomalous like-sign dimuon charge asymmetry A^^ 
reported by the DO collaboration Q. The updated result is A^i = (-0.787 ± 0.172 (stat) ± 0.093 (syst))%, based on 
9.0 fb"^ data ^. The result is 3.9a larger than the Standard Model (SM) prediction of (-0.024 ± 0.003)% ^. This 
asymmetry is comprised by the wrong-sign asymmetries a^f for iJ^j ^ mesons [2l, IJ], 

A'li = (0.594 ± 0.022)af; + (0.406 ± 0.022)a^;. (1) 

From direct measurements by B factories [^], af; = —(0.05 ± 0.56)% does not deviate from the SM prediction ^. 
Imposing these two experimental values into Eq. ([T]) , one finds a large a^j . The very recent update used muon impact 
parameter to directly extract [2| 

af; = -(0.12 ±0.52)%, a^; = -(1.81 ±1.06)%. (2) 

The result of a^; is much larger than the SM prediction of (1.9 ± 0.3) x 10~^ Q- The current world average of a^;, 
done before the very recent update 0], is 01 

a^, = -0.0115 ±0.0061, (3) 

which is still much larger than the SM prediction. This anomalous result has drawn intense theoretical attention, 
including model- independent analyses [5|-|9|, and explanations from specific new physics models [10l - ll7l |. 
The wrong-sign asymmetry a*; can be derived from mixing parameters [l| 

s AT 2|ri2,.| . , ,., 

o-si = T tan05 = — sm0,, (4) 

AttIs lAms 

where the AF^ and Arris are the width difference and mass difference of Bs — Bg system, (f>s is the CP violating phase, 
and F12.S is the absorptive off-diagonal element of mixing matrix (see Section II. A for more detail). Note that a*; is 
bounded by 2|Fi2,s|/Ams. The short-distance calculation in SM predicts Q, 

AF«,SM = (0.087 ± 0.021) ps'^, 
AF,^sm/F«,sm =(13.3 ±3.2)%, 

Am,,sM = (17.3 ± 2.6) ps-i, (5) 

(/.,,SM =(0.22° ±0.06°). 

Note that 0^ is very small in SM, so 2|Fi2.s| — |AFs|. If one inserts Eq. ([5]) into Eq. ([4]), one gets the small value 
of a^j mentioned before. These mixing parameters can be measured independently. In particular, Atti^ has already 
been well-measured. The current world average is [J] 

Arris = (17.78 ± 0.12) ps-\ (6) 

which is consistent with the SM prediction. Using the experimental Arus and a^;, Eq. ([4]) shows that Fja.s has to 
be enhanced by at least 3 times of the SM prediction. In fact, one of us has already pointed this out [18| in 2007, 
based on the earlier result of DO, which has almost the same central value as Ref. [l|] but with larger uncertainty. 
Recent studies fs, i?| also indicate this problem. On the other hand, AFg and (j)s can also be measured in several ways, 
although the precision is not as good as Arris- One method to extract these values is to study the Bg — >■ J/i^fj) decay. 
DO [3 reported 



AF, = +0.15 ± 0.06 (stat) ± 0.01 (syst) ps 
0s = -0.76+[5:3^ (stat) ± 0.02 (syst). 



+ 0.38 / , ,N , « rir, / ,^ '-'-' 



using 6.1 fb of data. The consistency of data between mixing parameters {Arris, AFs, and (fis) and a^; has been 
observed [J, l5| . Using almost the same amount of data, CDF [20| assumes (/)s = and reported 

AF, = +0.075 ± 0.035 (stat) ± 0.01 (syst) ps'^ (8) 

This central value drops to half the DO result, even below the SM prediction. But the two results still agree with 
each other because the uncertainties so far are still large. The consistency hints that new physics may play a role in 
Bs — Bs mixing. New physics can easily enter the dispersive M12.S and the phase <f>s- On the other hand, F12.S is 
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FIG. 1. The diagrams of Bs and Bs decay to Ds* Ds modes 



absorptive and thus hardly affected by new physics at high energy scale. As very many properties of B mesons have 
been studied and found to agree with SM predictions, new physics has to be rather exotic to change ri2^s while not 
affecting other known properties appreciably. 

The absorptive nature of ri2,s also makes the theoretical calculation challenging. It is helpful to revisit the calcula- 
tion of ri2,s in SM. One either approximates AF^ by operator product expansion (OPE) in short-distance picture, or 
estimates AFg from several modes which are believed to be important. The SM prediction [3] mentioned previously 
adopts the short-distance scheme. On the other hand, Aleksan et al. yj\ estimated AFj from exclusive two-body 

decays, mainly Dl Ds modes through color-allowed diagrams, as depicted in Fig. [1] Their result is close to the 
current SM prediction. They further pointed out that AFg induced by Ds ' D), modes approaches the result of 
parton model when the limits {mi, — 2mc) — >■ 0, rric — > oo and the large N^. limit are simultaneously imposed (for a 
detail discussion, see Ref. [221 )• How well does such an approximation hold in Nature remains to be checked. For 
example, as Ref. Q and one of us [l^ have already pointed out, a 100% long-distance correction is possible. The large 
all therefore motivates one to investigate the long-distance effect. In this paper, we perform a detail estimation of 
ATs from hadronic modes, which includes the two-body modes £)^^*-'**>jjy'**' ^ D^*'**^ Dsj{2700), and the three-body 
£)(*)£)W^(*) modes. 

We give the first estimation of the contribution to AFs by three-body Ds* D^*'>K'-*'> modes (and their CP conju- 
gates) . We use factorization approach, which seems to work well in color-allowed charmful three-body decays ^^ , in 
our calculations. As shown in Fig. [51 these modes can be produced by the diagram in Fig. [TJ but with an extra qq pair 
produced either in the current or in the spectator part, which we denote as current-produced {J') or transition (T) 

modes. The number of Ds'D^*^K channels are four times larger than Dg Dg modes, with a factor of two coming 
from extra qq^ which can be uu or dd, and another two from the choice of qq in either current or transition processes. 
With this enhancement in number of modes, and if the branching fractions of these modes are not very small com- 
pared with Ds Ds modes, it is natural to expect that AFg may receive non- negligible contributions from three-body 
jjy) Jj{*)xi*) modes. So far, the available measurements on these three-body modes are limited to current-produced 
modes with K in B^d systems only [24| - [28| . These modes are related to the corresponding modes in Bs system under 
SU(3) symmetry. We need to reproduce existing three-body data, before we make predictions for the Bs modes. 

Let us briefly survey the experimental situation regarding the SU(3)-related three-body modes. There is no measure- 
ment of either transition modes or modes with K*. Despite a 2.2a discrepancy on the branching fraction of B~ — > 
D'^D'^K~ decay between measurements [2^ [23], the branching fractions of current-produced Bu^d — >■ D^*'^D^*'>K 
modes are around 10~^ to 10~^, one order of magnitude smaller compared to two-body modes. So far, cs resonances 
i)gi(2536) and L'gj(2700) have been observed in the decays i?„,d -^ D^*^^D'^*^K (H-i^. For i:>gi(2536) resonance, 
its contribution to the branching fractions of three-body decays is in the order of lO"**, which is small compared 
with the total branching fraction. On the other hand. Belle observed that Ds,j{21Q{)) contributes to about half of 
the total branching fraction of B^ -^ D^D^K^. Note that L'sj(2700) has a fairly broad width (-- 0.1 GeV). These 



1 Throughout this work, we use D** to denote D*o(2317), Dsi(2460), or Z)si(2536). 
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FIG. 2. The current and transition diagrams. The left part, Bs -^ Dl* D^*' K^*' , is the current-produced diagram, and the 
right part, Bs — > Di* D^*' K^*' , is transition diagram. 



measurements suggest that £'si(2536) could be treated in a two-body picture while it is more appropriate to con- 
sider Dsj{2700) in three-body decays. Furthermore, the contribution of 1)^,7 (2700) in B~ — > D^D'^K^ decay is 
B{B- -^ D°D,j{2700)) x B{D,j{2700) -^ DK) = (0.113to:o4o)%, which is about half the total branching fraction 
(0.222 ± 0.033)% '26f|. Consequently, the contribution of Z),,j (2700) in three-body modes and in AF^ should be 
investigated. 

This paper is organized as follows. In Section |TT1 we describe our formalism and briefly review the newly discovered 
Z?s J (2700) resonance that has a non-negligible contribution to three-body modes. The results of two-body modes are 
in Sec, nil Al For three-body modes, we examine the factorization formalism and calculate AFs in Sec. lIIIBl Another 
scenario and the effect of four-body modes are discussed in Sec. IV, followed by the concluding section. Numerical 
inputs and some calculational details are collected in three Appendices. 



II. FORMALISM 



A. Formula for AF 



The time evolution of a B^ meson can be described by the following formula. 



= M 





(9) 



in which we adopt the phase convention of \B) and \B) to be CP\B) = — |i?).G The F term in Eq. © is the absorptive 
part, which can be calculated by summing all on-shell intermediate states. 



2Afe- 



j d^A%^f{^)AB,^fm, 



(10) 



where $ is over phase space variables. 

We define the width difference AF^ as the difference between light and heavy eigenstates, F^ — Th- Assuming CP 
conservation, which is a good approximation for SM in Bg — Bs system, the eigenstates of Bg meson are CP even and 
odd states. From short-distance calculation of SM, the light and heavy eigenstates correspond to CP even and odd 



Our phase convention differs from that in Ref. [i||. 



^ For n-particle mode, the phase space measure is d<l? = 0?=! Tg^ ^ (27r)''<5*( \J pj 
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states respectively. Thus, the AF^ can be related to Fy by 



AF EE Fl - Th 
= - 2Fi2 



B ^ J 



(11) 



in which we have used F21 — F*2 from CPT symmetry, and F12 = F^2 = Rc(Fi2) from CP symmetry. The fact that 
F12 is real under CP symmetry can be seen from 

B j: J 
B j: J 

1-5] Re J d^A*s^f{<^)As^fm . 
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The amplitude product Ag^f{^)A§^f{^) is complex conjugate to the amplitude product of conjugate intermediate 
state A*g^A^)AB^f{^) by CP symmetry. F12 sums up all the intermediate states and turns out to be real. For 
convenience, we define the width difference of each exclusive decay as AF/ , and its corresponding complex term in 
F12 to be 



AF 



/ 



where F12 j is defined as 



12,/ 
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-2 X Re[Fi2,/] 



d'^A*B^fi'^)As^fm. 



(13) 



(14) 



Although F12./ is complex by looking at one mode, the imaginary part is cancelled by its CP conjugate mode, and 
thus the total F12 turns out to be real. Once As-^fi^) and Ag^j:{^) are known, one can readily calculate the 
corresponding AF/ and branching fractions. In the next section, we will apply the factorization formalism to obtain 
these amplitudes. 

Before we move to model-dependent calculation, it is useful to extract some general limits of the magnitude of AF/ 
from Eq. (|13l) . For an intermediate state |/), the magnitude of AF/ induced by this state is bounded by 



(15a) 
(15b) 
(15c) 
(15d) 



AF/ 
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= f ■ |Re[2Fi2, 


<^-|2ri2./| 
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^ / d^,/\A~^.J\A^ZA^) 
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B^f- 



There are three inequalities in this formula. The first inequality reflects that AF/ is only proportional to the real part 
of ri2,/- The second inequality is obtained by the fact that the phase of the amplitude product A*g^f{^)AB^f{^) 
may be different over the phase space, which would reduce the overall |Fi2,/|. For the last inequality, it accounts for 
the "mismatch" effect between |^b^/($)| and |y^3_^/(<i>)|. Even when the branching fractions of i? — >■ / and B — >■ / 
are the same, the induced AF could be quite small if the decay probabilities of the two modes are highly mismatched 
in phase space. Note that the latter two limits are experimental observables. If the branching fractions oi B ^ f and 
B ^>- f are measured, one could find the maximal magnitude of the corresponding Fi2j. The bound can be refined 



by the second inequality if the Dalitz plots of the two modes are available, 
real part F12,/, could be any value in the range of — 2|Fi2,/| to -|-|2Fi2j|. 



But the AF/ , which is proportional to the 



B. Factorization Formalism 

The relevant effective Hamiltonian for the b -^ c transition is 

n.s = ^yefcV;:[ci(A*)OJ(A^) + C2{fi)0^,if,)] , (16) 

where Ci{fj,) are the Wilson coefficients, and Vcb and Vcs are the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments. The four-quark operators Oi are products of two V — A currents, i.e. Ol — {cb)v^A {scjv-A and O2 = 
{sb)v-Aicc)v~A- 
With the factorization ansatz, the ampUtudes for two-body Bg — > Vl Vl decays are given by 

A{Bs ^ I?(*)P(*)) = ^V,tV;^a,{Vi*M{cb)v-A\Bs){'Di*'>\{sc)v^A\0}, (17) 

where the effective coefficients are expressed as ai = ci + cijZ if naive factorization is used. Note that Vs could be 
the usual D^*^ and or higher D^ resonance such as £»*o(2317), i:»^i(2460), i:»si(2536), and Dsj(2700). The factorized 
amplitudes consist of the products of two common matrix elements: the current-produced uk and the B^ to I?s 
transition. They are parametrized by the standard way [2^. The matrix elements of current-produced X>s are 

(18) 



(I?,(p)|(y-A)jO)= */p,p^, 
(2?:(p, A)|(F - A) JO) = mc./i,.e;(p, A). 



The transition matrix elements for T)s are 



,,„^, 2F3^°^-(g^) 



{V:{pD^,X)\{V-A)^\B,ipB)) = e^,p.£*»^. • 



ms + mo* 



g2 



' '%J {mB + mD.)F^'^'{q^) (19) 



/ Tfl TTl \ 






-* — :^qf,2mD'FQ ' (q ), 
q'^ 

where e^j^po- is the totally anti-symmetric symbol with £0123 = 1- For convenience, our notations of decay constants 
and form factors of D** are different from the usual notations. The conversion can be found in Appendix IB] 
The amplitudes of three-body modes £)(*)z)(*)/^(*) decayed from B and B are given by 

Aj{B, ^ i?(*)5(*)i?(*)) = ^V,bV:MDi*Mcb)v-A\Bs) ■ (5(*)i?(*)|(sc)y_^|0), 

(20) 
At{Bs ^ i?i*)5(*)i?(*)) = ^KbK>i(5(*)i?(*)|(c5)v_^|i?,) • (Z?(*)|(sc)y_A|0), 

where y^j- and Aj- denote the amplitudes of current and transition diagrams, respectively. Unlike the Vl 'Vl ' modes 
in which only standard form factors appear, these amplitudes involve the time-like form factors and transition form 
factors of two pseudoscalars (DK) or vectors {D*K*), or a pseudoscalar with a vector {D*K or DK*). 

The parametrization of time-like form factors are similar to the space-like counterparts, such as {Ds* \V — A\Bs). 
The time-like form factors of two pseudoscalars (PP) states are given by 

{PaiPa)PbiPb)\iV - AUO) = (^ipa-p,), - ^^^^^Q,) Ff^ iq') + '^^^q.Friei (21) 



where q^ — Pa + Pb ^^ ^^^ momentum of the current. For the states with one vector and pseudoscalar {VP), the 
parametrization of time-hke form factors are 

{V{pv,ev)P{pp)\{V - A)^|0) = - 6^.p.£t.»^ • '^^^^ - * (^y^ " ^«^) ('"^ + nip)AX''{q^) 



« ( [Pv-Pp),. -^ — -q^i (ey •'7) 7 i^-^^q^,2mvA^'^ {q^). 



(22) 



The time-hke form factors of two vectors {VV) states can be parameterized analogously, 



[■nia + mbY [ma + mb) 



{Va{pa,ea)Vb{pb,eb)\{V - A)f,\Q) = if-aypcrel"'el''p''^plqt, iJ' ^ _ + if-t,ypas7p''aPl{£l ■ q)-, x^ 



+ ie,.,.eiy^^pl{sl ■ ci)^^^, + (<, - ^g.) {el ■ q)AViq') 



{^l, - ^?m) K ■ 9)^r2^(?^) + (iPa Pb)p "^^^^Ip) {el ■ eDAX'^iq') 

(23) 



+ {el ■ q){el ■ qf^AXr{q% +{el ■ el)^{ma + mbYA^.'^iq') 



The transition form factors are more complicated. The case of B^ to PP transition form factors were formulated 
in a general way in Ref. [301] , which can be rewritten as 

{Pa{Pa)Pb{Pb)\{V ~ A)^\Bs{pb)) = e^.paPW{Pa - PbT ^ + i [{PB + q) p ^ ""^^ g^ ^ 

+ * (Pa - Pb)p -^^<lp — + ^ 2 gA' 



q J ms, g' mB, 



where q'^ — p^ + Pb is the total momentum of PP, and g'^ = p^ — g'^ is the momentum of the external current. In 
this form, the terms with Ai and A2 are zeros when contracted with q' and q. For the transition form factors of Bg 
to VP and VV, since they are more complicated and there is so far no data, we only write down the form factors 
obtained from pole model rather than the general forms. For VP, we have 
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FIG. 3. Pole diagram of Bs ~ Bs mixing. Left: current-produced diagram. Right: transition diagram. 
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and for Bg to VV , we parameterize 
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Under CP conservation, all 
transfornrations are provided i 



ml -q^ ,\ Af/^^ g' A^f^^ 
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(26) 



P conservation, all these form factors can be related to the form factors of their CP conjugates. These 
ions are provided in Appendix |B] 



C. Pole Model 



Since the branching fractions of Dg Dl are large, it is natural to expect a sizable contribution from off-shell 
Ds poles. In addition, experiments have observed -D^j (2700) in the three-body decays as we have described in the 
introduction [23,[2g|. Dsj{27QQ) can decay to on-shell D^*''K, but only goes off-shell to D^*'^K* because of kinematics. 
As shown in Fig. |3l we consider pole exchanges, including Dg, D* and Dsj{2700), in three-body decays. Note that 
the Ds pole only goes to D*K rather than DK. 

In the following calculation, we use off-shell Ds' poles and -Dsj(2700) to model the D'-*'>K^*^ form factors. The 
effective Lagrangian taken from Ref. |3ll - l33| is applied to describe the interaction between Dq*' mesons and light 



psedudoscalar or vector mesons. The pole contribution to D^*'^K^*^ form factors can be calculated by 

(7?Wi^W|(l/-A)JO)polc=^ ^^ X (i^Wi^WK/:eff|Ant)(A„t|(l^-A)^|C 



9^ - "^^nt* + *"T-int*rint, V ' "^Lt* 






;i?wxW|z/:efflA„t>(A;tl(^-^)Jo) 



9 ^ "^mt + «"iintri„t 



(27) 






1 ""int* 



'D(*^^K^*'>\^CMD*^,){DU{V - A),\B)) , 



where the D\J^ is the intermediate particle with mass JTiintf*) f^nd width Tint^^). We adopt the Breit-Wigner form of 

the propagator and replace efj^ej^^ as (— g"'^ + g"'Z'^/'Tif„^^) to account for the off-shell effect. The explicit forms of the 

matrix elements {D'^*'^ K'^*'^\iCcs\D-*-l) in the above equations can be found in Ref. [33|. A full list of pole contribution 
to form factors are listed in Appendix [Cl 



D. Dsj(2700) Resonance 

The relevant properties and parameters of 13^,7 (2700) are summarized in this section. The mass and width of this 
resonance are [34| 

mD.,(2700) = 27091^ MeV, 
rD,.(27oo)-125±30MeV. 

Note that the width has a large uncertainty (^ 25%). The ratio of branching fractions of this resonance to DK and 
D*K \s also measured 



^ g(A.(2700)+^ D*K) ^ 

where D'^*^K is the average of D^*^Ks and D''*^K^ modes. On the other hand, the contribution of Z?sj(2700) in the 
decay B+ -^ D^D"K+, denoted as B{B+ -^ D°Dsji2700)) x 6(A,/(2700) -^ D°K+), is extracted [2g| 

B{B+ -^ D^Dsj{2700)) x S(Aj(2700) ^ DK) = (11.3l^;|^) x 10"^ (30) 

which constitutes about half the total branching fraction of this measurement. Note that this quantity has a large 
uncertainty, similar to the measurement of width. The quantum number of I?sj(2700) is determined to be J^ = 1^ 
from helicity angle distribution, which limits this resonance to be either an s-wave or d-wave meson (or a mixed state 
between them). The interpretation of I?s,/(2700) as a radial excitation of D* [n^^^^Lj = 2^Si) is proposed, which 
can explain its mass 36[, partial width [37], and contribution in 5+ — >■ D^D^K'^ decay [38|. In some strong decay 



models, a mixed state 2'^ 5*1 — l^Di describes the partial width better [39|. As the theoretical predictions of mass and 
partial width are highly mo del- dependent, the identification is still not clear yet. We assume Z3sj(2700) as a 2^ Si 
state in this study. 

The effective Lagrangian in Ref. |3ll-[33l| can still be applied to describe the interaction between Dsj (2700) and light 
mesons [33|. We work out the relevant matrix elements, 

{D{p2)K{p3)\i£cs\Dsj{2700){pi,ei)) =-igD,.,DK£i ■ Ps, 

a I O J- 1 

{D*{p2,£2)K{p'i)\iCcs\Dsj{21QQ){pi,£i)) ^ - igo^jD- K et,ual3£i£2P3Pl: 
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Mode(f) 


D°K+ 


D+K'' 


D*°K+ 


D*+K° 


DsTI 


D'sT^ 


r{f) 


1.02 


0.98 


0.93 


0.89 


0.17 


0.04 



TABLE I. The ratio r of the branching fractions of six main decay modes of the Dsj(2700)+ resonance. 



where the strong coupUng constants are given by Ref. i37i] . 



^: ^ (32) 

o .9 jmo' 
gOsjD'K — ^-f- 



h V ^D,j ' 

with /tt = 132 MeV. Once the couphng constants and form factors are extracted, one can insert Eq. (pij) into Eq. (P7)) 
to obtained the contribution to form factors from the Dsj (2700) resonance. 

From these matrix elements, Ref. [37] predicted the ratio of branching fractions 

r(i:)*AO = 0.91 ±0.04. (33) 

This ratio agrees with Eq. (j29p very welL The ratios of the branching fractions of the six main decay modes are given 
by Tabic HI The mixing angle between 77 and 77' is taken from Ref. [40| . 

Assuming £'sj(2700) only decays to D^*'>K and D'-*^rj^'\ g^ is proportional to the total width. Thus, we have 

5 = 0.28 ±0.03, (34) 

where the uncertainty comes from the uncertainty of the total width. Note that this value is slightly larger than the 
one in Ref. [33| as the world-average of width [Eq. ([25)1 ] became larger. 

Taking the measured mass, width and B{B -^ D''*^Dsj {2700)) x B{Dsj{2700) -^ DK) (see Sec. Ill B for details) 
as input, the Z?s,/ (2710) decay constant is extracted as 

/D.,(27oo)=240±31MeV. (35) 

The decay constant can be compared to the previous estimations 243 ± 41 MeV in Ref. [33| and 295 ± 13 MeV in 
Ref. 133|. Note that it is compatible to the decay constants of Ds\ which we use 260 ± 13 MeV in later calculation. 
The Bs — > i'sj (2700) transition form factors can be obtained by using a covariant light- front quark model |33j . 
For the 25 wave function, □ its Gaussian width can be fixed by the decay constant derived from Eq. ([50)1 . It is then 
straightforward to obtain various Bg — ^ Dsj form factors: 



^B.Z5.,(2700)(^2) 
^B.Z5.,(2700)(^2) 



0.25 ±0.03 



1 - 0.03 gVm|^ + 0.38 q^/m%^ ' 
0.24 ±0.02 



0.17 ±0.02 ^ ' 



1 + 0.66 gVml^ + 0.54 q^/m%^ ' 
0.007 ±0.001 



.BsD,j(21W), 2-. ^ 

These transition form factors are small comparing to the Ds (collected in Appendix [A| . because of the poor overlap 
between wave functions of ground state B mesons and the radial excited Dsj(2700). 



** In the quark model with a simple harmonic like potential, the wave function for a state with the quantum numbers (n,l,m) is given 
by f„i{p^//3'^)Yijyi(j>)ex.p{—jP/2l3) with /io(x) = 1 and f2oix) = V 4(~1 + 3^)' ^^ ^^ ^^'^ Gaussian width /3 to decay constant. 
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Mode(f) B(B,,(„) ^ /) (%) B{Bs ^ /) (%) B(B.^ /) (%) Ar^/r^ (%) AF^/F^ (%) 

data this work Ref. [21j| this work Ref. [23| 



Ds-D« 1.04 ±0.35 

(1.00 ±0.17)0 

DlDs+DsDl 2.75 ±1.080 
(1.58 ±0.33)0 

DlDl 3.08 ±1.490 
(1.71 ±0.24)0 



1.4 ±0.3 ±0.3 1.6 



1.8 ±0.4 ±0.4 2.2 



2.3 ±0.5 ±0.5 3.6 



2.7 ±0.6 ±0.6 3.1 



3.6 ±0.8 ±0.8 4.4 



3.8 ±0.8 ±0.8 6.9 



Di*'_Di*^ 4.9 ±1.4 



6.9 ±2.3 
4.0 ±1.5 
(4.29 ±0.74)0 



5.5 ±1.2 ±1.1 7.4 



10.2 ±2.2 ±2.1 14.4 



TABLE II. The branching fractions of Ba -^ Ds Di* decays and their contribution to the width difference. The results can 
be compared with data in Refs. [41. |34|. l4l| . The data for B~ system in Ref. I34[]. which are related to Bs under SU(3) symmetry, 
are shown in parentheses (see text for detail). The theoretical result of Ref. [2l| ] is also presented for comparison. 



^ Data taken from Ref. 
^ Data taken from Ref. 
'^ Data taken from Ref. 



E. Non-Resonance Contribution 

In general, there will be both resonant and non-resonant (NR) contributions to form factors. In previous study 
of B ^>- D^*'> K^ K^ decays [23], it is necessary to add NR contribution to form factors to explain the experimental 
observations. Therefore, we should include the NR effect in this work. To produce the d'-*^K^*'> pairs, at least 
one gluon must be emitted to produce qq pairs. The QCD counting rule [23| provides an ansatz for the asymptotic 
behavior of the non-resonant form factors, which is 



F{q^)NR 






In 



A2 



(37) 



where q^ is the invariant mass of d''*^K'^*^ and A = 0.5 GeV is the QCD scale. 

Together with the pole contribution provided in Appendix [Cl the complete form factors are modeled by the pole 
and NR contribution, 



F(g2) = F(Q2)|p,,, + ^ 



In 



A2 



(38) 



where the asymptotic form of NR contribution is adopted for simplicity. As more data is available in the future, one 
could replace this simple form with a more sophisticated one to fit the data, as in Ref. |23j . 



III. RESULTS 



A. Two-body I?i*'oi*' Decays and the Width Difference: An Update 



We first update the branching fractions of two-body Bg — > Dl ' Dl ' decays, which contribute to AF^. The 
necessary parameters are given in Ap pen dix \K\ Our results are listed in Table [Til where experimental results and 
previous theoretical results from Ref. [21[ are listed for comparison. Since SU(3)-related modes in Bu,d systems are 
usually more precisely known than the B^ system, we also list them in parentheses for comparison. For example, data 
for B{Bu — ^ DuDs), which is approximately the same as B{Bs -^ DgDs) under SU(3) limit, is listed. Note that two 
uncertainties are given in our results: The first uncertainty is obtained by varying decay constants and form factors 
by 5%, while the second comes from the estimated 10% uncertainty in ai. 

The branching fractions of Dg Dg modes are all of percent level. In general, our result is smaller than the result 
in Ref. [2l|. These branching fractions can be compared with experimental data in both Bg and B~ system. One can 
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Mode(f) B{B, -^ f) (%) BjBs -^ f) (%) AFj/r, (%) 

D._D:o(2317) 0.10 ±0.02 ±0.02 0.15 ± 0.03 ± 0.03 -0.24 ± 0.05 ± 0.05 

(0.073lH??)Q 

D* 0*0(2317) 0.05 ±0.01 ±0.01 0.12 ±0.03 ±0.03 -0.15 ± 0.03 ± 0.03 

(0.09 ±0.07)0 



D,Dsi{2460) 0.24 ±0.05 ±0.05 0.04 ± 0.01 ± 0.01 -0.18 ± 0.04 ± 0.04 

(o.3il«:S^) 

L»: 5^1(2460) 0.81 ±0.17 ±0.17 0.06 ± 0.01 ± 0.01 ±0.16 ± 0.03 ± 0.03 

(1.20 ±0.30) 

D.L>,i(2536) 0.02 ±0.01 ±0.01 0.38 ± 0.08 ± 0.08 ±0.19 ± 0.04 ± 0.04 

(0.022 ±0.007)0 

D*5,i(2536) 0.09 ±0.02 ±0.02 0.38 ± 0.08 ± 0.08 ±0.34 ± 0.07 ± 0.07 

(0.055 ±0.0016)0 

£):o (2317)5^1 (2460) 0.024 ± 0.005 ± 0.005 0.002 ± 0.001 ± 0.001 ±0.013 ±0.003 ±0.003 

7?:o(2317)Osi(2536) 0.002 ± 0.001 ± 0.001 0.017 ± 0.004 ± 0.004 -0.012 ±0.003 ±0.003 

73rf (2460)5^1 (2536) 0.001 ± 0.001 ± 0.001 0.077 ± 0.017 ± 0.016 ±0.000 ± 0.000 ± 0.000 

i3:o(2317)_D:o(2317) 0.009 ± 0.002 ± 0.002 ±0.018 ± 0.004 ± 0.004 

i3si(2460)5»i(2460) 0.014 ± 0.003 ± 0.003 -0.010 ±0.002 ±0.002 

73rf (2536)5^1 (2536) 0.007 ± 0.002 ± 0.001 ±0.008 ± 0.002 ±0.002 

Total 2.57 ± 0.55 ± 0.54 0.24 ± 0.27 ± 0.05 

TABLE in. The branching fractions and width difference of Bg and Bg decays to two-body D^* '** I)** , where D** is Z)*o(2317), 
Dsi(2460), or Z)si(2536). We show data of SU(3) related modes in Bu system [3J| in parentheses for comparisons. 

^ B(B~ -> D(*)OSso(2317)) X Z3(Sso(2317) -^ Dsir-). 

^ B{B- -> D(*)0_Dsi(2536)) x B(_Dsi(2536) -> p*K-). 

'^ The contribution from CP conjugate modes / is included. 

see that our results agree with experiment within uncertainties. The direct measurement of Bg — >■ Dl Dl exclusive 
decays was recently reported by Belle [41j. While the observed branching fraction of DgDs mode (1.0 ± 0.4)% is 
close to our result, other modes are more aligned with the calculation in Ref. 21]. But the world average of the 

inclusive branching fraction B{Bs -^ Dl* Ds') [J, [SJ and the rates of SU(3) related modes are closer to our results. 
The total ATj/Ts induced by Ds* Ds* modes is 10.2 ± 2.2 ± 2.1%. This value is smaller than the previous long- 
distance calculation [2l| also shown in this table. In addition, the total ATf/Ts does not reach the short-distance 

central value in Eq. ([5]). One also observes that ATsiDs Ds )/rs is approximately two times the total branching 
fractions. The relation \Ars{f)/Ts\ < 2yjB{Bs -^ f)B[Bs — >• /), which corresponds to the maxima in Eq. (J15d|) . 
saturates only when the mode(s) / is purely CP-even, such as the DgDs mode. The nearly maximal ATf reflects 

that Ds are very efhcient in mediating the width difference. 

Several new cs resonances are found in B decays. They may also contribute to AT^. We calculate the contribution 
by the two-body modes with L»*o(2317), L'si(2460), and i:'si(2536). Results are shown in Table |llll There are 
additional 21 modes when these higher D** resonances are included. Note that not all modes are shown explicitly in 
the Table. Since CP is conserved in this work, B{Bs —>■/) = B{Bs — >■ /) and ATf = AFj. For modes which are not 
CP eigenstates, the contributions from their CP conjugates are also known and should be added to ATs/Ts. The total 

branching fraction of these additional modes is comparable to the sum of B{Ds Ds ). However, the corresponding 
contribution to the width difference turns out to be tiny. After considering all of these two-body modes, the total 
ATf/Ts only increase slightly from 10.2 ± 2.2 ± 2.1% to 10.4 ± 2.5 ± 2.2%. There are two reasons for such a tiny 
contribution. First, the sign of ATf are fluctuating among these modes, leading to cancellations in the total sum. In 
addition, the "mismatch" effect is serious. For instance, the Bg — > -D*5si(2460) mode has a non-negligible branching 



^ Note that this measurement does not tag the flavor of the Bs meson. Although there should be a corresponding correction to the order 
of ATs/Ts [22l| . it is smaller than the theoretical errors and omitted from the table. 
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Mode(f) BjBs -^ /) (%) BjBs -^ f) (%) AFj/r, (%) 

DsDs.,{27m) 0.44 ±0.18 ±0.09 0.02 ± 0.01 ± 0.01 0.21 ± 0.08 ± 0.04 

L»:L',j(2700) 2.0 ±0.8 ±0.4 0.08 ± 0.03 ± 0.02 0.73 ± 0.27 ± 0.15 

Di*'L',,7(2700) 2.5±1.0±0.5 0.11 ± 0.03 ± 0.02 1.9 ±0.7 ±0.40 

DTDsj{270Q) 0.14 ±0.08 ±0.03 0.02 ± 0.07 ± 0.01 0.08 ± 0.03 ± 0.02 

TABLE IV. The branching fractions and width difference of two-body Bs and Bs decays to Da*'** Daj{27QQ), where D** stands 
for D:o(2317), Dsi(2460), or _Dai(2536). 

" The contribution from CP conjugate modes / is included. 

fraction 0.81%, but the branching fraction of B^ — > D*D si{2AQQ) is only 0.06%. In fact, the smallness of contributions 
in the heavy quark hmit from p-wave resonances was expected [2l| , and is confirmed in a reahstic calculation given 
here. 

The sizable branching fraction B{B -)■ D^*'>Psj (2700)) x B(Z)s,/(2700) -^ DK) indicates that the D,j{2700) 
resonance may be important for AF^. Since Z?sj (2700) has a broad width, it is expected to interfere with the 
continuum of Bs — >■ DsD^*^K produced by Dg* poles (see Fig. |3] and the next subsection). For completeness, it is 
better to calculate the contribution of Dsj{2700) to AF^ in three-body modes, including the on-shell and off-shell 
parts. However, the two-body calculation is simple and straightforward. It is, therefore, helpful to see the contribution 
of Di*'**'^ Dsji2700) to AF^ first. 

Using the parameters calculated in Eq. (p6|. the contributions from two-body modes including Dsj {2700) is shown 
in Table ITVl Several things ought to be noted: (a) The branching fractions of modes with current- produced Dgj {2700) 

(the B{Bs — >■ /) column of Table. IIV[) are comparable to those of the Dg Dg* modes. The two-body decays with 
Dsj{2700) seem to be suppressed seriously by phase space at first glance. Nevertheless, this may not be true since 
the factorized amplitude {'Dl\{V — ^)^|0) [see Eq. (fTSl) ] for current-produced meson is enhanced by mass, and the 
decay constant of Dsj{2700) is unsuppressed. (b) For the mode Bs — > D*Dsj{2700), there are several enhancement 
and suppression factors, when replacing D* with Dsj {2700). First, its amplitude is dominated by s-wave and is 
free from additional momentum suppression. In addition, it is enhanced through the above mentioned factorized 
amplitude and suppressed by phase space. The branching fraction of Bg -> D*Ds,j {2700) turns out to decrease 
~ 10% compared with Bg -^ DID*. On the contrary, the decay Bs — >■ -0^-0^,7(2700) is p-wave. Its amplitude and 
thus branching fraction drops more than 50% when compared to Bg — >■ DgD*. The two different trends lead to a large 
ratio B{Bs -^ D*s Ds.j {2700)) /B{Bs -^ DsDsj{2700)) w 5. (c) The branching fractions of modes in which Dsj{2700) 
contains the spectator quark (the B{Bs — > /) column) are very small. The branching fractions are suppressed not 
only by phase space, but also by the small transition form factors shown in Eq. p6p . 

_ The ATs from Di*^Dsj{2700) is 1.9 ± 0.7 ± 0.4%. As the upper bound in Eq. (flSdl) implies, the AF JF^ of 
-Dsj(2700) is limited by the imbalance between the modes in which -Dsj(2700) produced via current or with spectator. 
Nevertheless, the contribution form -D^j (2700) is larger than those from D** and should not be neglected. We remark 

that, as we shall see in the three-body case, the transition amplitudes from Dg poles can interfere constructively 
with the current-produced Dgj pole and overcome the above mentioned suppression, leading to sizable contribution 

to AF,. 



B. Three-body oi*^ D'-*'> K^*^ Decays and Contributions to the Width Difference 

We now turn to the three-body case. We shall first compare our results with the measured branching fractions in 
Bu,d system, starting from pole model and including NR effect, if necessary. After demonstrating that our calculation 
is consistent with data, we proceed to calculate the width difference in the Bs system. 

1. Current-Produced Branching Fractions in B^.d systems 

Only current-produced modes with K have been measured in B^^d systems. There is no measurement for the 
rest of the modes, including current-produced K*, and all the transition modes. A summary of current data and 
our results is presented in Table |Vl We separate the results of BaBar and Belle for comparison. Note that in B^^d 
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Measurement 



BaBar Data(%) Belle Data(%) Our Results (%) 

Scenario I (I') 
Pole model with Dsj 
(without Dsj) 



Remarks 



Scenario II 
Pole model+NR 



Category 1: current-produced DK with B ^ D transition 
B{Bu-^DuD,j{2700)-)x N/A 0.113J:J!:|54^ 0.12 ± 0.08 ± 0.03 0.12 ± 0.08 ± 0.03 Input for Scenario I 



BiDsj{2700)- ^D°K' 



(0) 

0.131 ±0.0143 0.222 ±0.03^ ~ 0.23 

(~ 0.07) 



B{Bu -> DuD"K-) 

B{Ba^DdD°K-) 0.107 iO.OllQ N/A 

B(Bi^DiD,.,{27m)-)x N/A N/A 

B(D,j(2700)~ ^D°A-) 



0.11 



Color-suppressed di- 
agram neglected 



0.22 ± 0.14 ± 0.05 O.lOlJJo^ ± 0.02 Input for Scenario II 

(0.06 ±0.03 ±0.01) 

0.11 ±0.07 ±0.02 0.11 ±0.07 ±0.02 

(0) 



Category 2: current-produced DK with B -^ D* transition 

0.247 ± 0.0210 N/A 0.67 ± 0.45 ± 0.14 0.32tol3 ± 0.07 Input for Scenario II 

(0.07 ±0.03 ±0.01) 



B{Ba-^Dl&'K-) 



B{Bd -^ D'aDs.j{2700)-)x N/A 
B(D,j(2700)- -^D°K-) 



N/A 



0.50 ±0.33 ±0.11 0.50 ±0.33 ±0.11 
(0) 



BiBa~^DaD*°K-) 



B{Bd -> DdDs.j{2700)-)x N/A 
B(D,j(2700)- ^D*°K-) 



Category 3: current-produced D* K with B ^ D transition 

0.346 ± 0.0410 N/A 0.35 ± 0.21 ± 0.07 0.35 ± 0.21 ± O.ofI 

(0.20 ±0.10 ±0.04) 

N/A 0.11 ±0.07 ±0.02 0.11 ± 0.07 ±0.oS 

(0) 



BiBa -^ D'aD*°K~ 



B{Ba -^ DID S.J {27m)- )x N/A 
B(D,j(2700)~ ^D''°K-) 

B{Ba -^ DlD*+K'') 



B{B^ ^ D'^D*+K%) 



Category 4: current-produced D* K with B ^>- D* transition 
1.060 ± 0.09^ N/A 

N/A 

0.826 ±0.08(0 N/A 

0.44 ± 0.0^ 0.34 ± 0.0^ 



0.94 ± 0.62 ± 0.20 
(0.15 ±0.08 ±0.03) 


0.94 ±0.62 ±0.20 






0.52 ±0.33 ±0.11 

(0) 


0.52 ± 0.33 ±0.] 


LlEl 






~0.91 
(~0.15) 


~0.9]I!| 




Color-suppressed 
agram neglected 


di 


~0.46 
(~ 0.07) 


~0.4^ 




Color-suppressed 
agram neglected 


di 



TABLE V. Comparison between experimental results from BaBar and Belle collaborations and our results in Scenario I, II, 
and I'. See text for detailed definition. 

^ Ref. 
^ Ref. 
" Ref. 
d Ref. 
'^ In Scenario II, the results of modes in Category 3,4 are the same as Scenario I. 



systems, some £)(*)z)(*)^(*) modes contain both color allowed and color-suppressed diagrams, where the latter is 
expected to be sub-leading and is neglected in this work. We labelled these modes in the remarks of the table, and 
also add approximation sign in front of our results. Note that in the calculation of AF^ in Bg system, color-suppressed 
diagrams only appear in modes with r]^'^ and do not affect D s D'^*^ R'^*^ modes. 

According to whether D or D* , there are four types of D'^*^D'--*^K modes, which are classified into four categories 
as shown in Table |Vl Modes in each category have similar branching fractions because of SU(2) symmetry. The 
measured branching fractions increase from Category 1 (~ 0.1%) to Category 4 (~ 1%). One can find tension in 
measurements of i?„ — ?► D^D^K^ . A large Z)sj(2700) contribution has been observed in B„ -> DuD^K~ by Belle 
only [2g|, but in 2.2cr disagreement with BaBar [28j. The tension in data becomes more severe if one compares the 
Dsj{27QQ) contribution to the total branching fraction of i?„ -^ DuD^K~ . In the case of Belle, the contribution 
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from _Dsj (2700) is about half the total branching fraction. However, it is approximately equal to the total branching 
fraction for BaBar. As we show, the inconsistency makes it difhcult to explain all data with a simple pole model. 

The results of our calculation in different scenarios are compared with experiments in Table |VJ In Scenario I, Dl 
and Dsj poles are used, while in Scenario r, only D), poles are considered, with results shown in parentheses for 
comparison. In Scenario II, NR contributions in DK time-like form factors are included to demonstrate that the 
inconsistency with experiments in Scenario I can be resolved. Note that no NR contribution is introduced for modes 
in Category 3 and 4 as the pole model results (Scenario I) already agree with data. Furthermore, as there is no 
measurements on transition modes and modes with K* , no NR contribution is applied to these modes. The two 
uncertainties of our results are obtained by the same method as in two-body case, but with additional uncertainties 
from strong couplings included in the first errors. 

Despite the disagreement between data, we first attempt to explain all measurements only with a pole model 
(Scenario I). The corresponding diagrams can be found in the left portion of Fig. [3] with the appropriate spectator 
quark. In the calculation, we first fix the decay constant of Dsj(2700) from the contribution of Dsj(2700) in iJ^ — >■ 
DuD^K~ decay. The value of this decay constant is shown earlier in Eq. p6p . and the value agrees with those 
obtained in other studies (see Section II. D). The total branching fraction of i3„ — >■ DuD^K~ is consistent with 
Belle's measurement, and inevitably less consistent with the BaBar result and the SU(2)-related mode Bd — >■ DdD^K~ . 
Unfortunately, there is no measurement on Bd — >■ DdD^K~ rate from Belle yet. For Category 2, the total branching 
fraction Bd — >■ D*^D'^K~ is about 2.5 times larger than the BaBar result as in Category 1. Again, there is no 
measurement from Belle. More data analysis is called for. Nevertheless, it is interesting to see that our predicted 
results on branching fractions in Categories 3 and 4 agree well with data. 

To explain the total branching fractions in Scenario I, we must start from the Z?sj (2700) contribution, which has 
on-shell as well as off-shell parts. Roughly speaking the i5sj(2700) contribution can be understood by using the 
narrow width approximation. The contribution in Category 1 (2) is almost the same as in Category 3 (4). This is 
expected since the two categories are different from each other only in Dsj{21QQ) -^ D*K, DK parts, which have 
nearly the same branching fractions [see Eq. (|33p ]. The contribution in Category 2 is about five times larger than in 
Category 1, where the B -^ D* transition is replaced with B -^ D. This factor already appeared in the two-body 
branching fractions of Bg — >■ Dg Dgj modes shown in Table HVl However, a closer look reveals that the precise 
Dsj{2700) contribution should be obtained by integrating the full three-body phase space, as the width of Dgj {2700) 
is of the order of 0.1 GeV, which is not narrow enough compared with the three-body phase space. (For instance, 
the decay Bs -^ D*Z)sj(2700) with Dsj{2700) -^ D*K, the invariant mass of D*K ranges roughly from 2.5 GeV 
to 3.3 GeV. The Breit-Wigner function for Dsj{27QQ), with a peak at 2.7 GeV, cannot be approximated as a delta 
function since its peak is less than 2 times of width above the lower limit of the invariant mass of D^*'K.) The 
numerical results usually show a 10% overestimation by narrow width approximation. In addition, the Dsj{2700) 
contribution in Bd -^ D*j^D*^K~ is shghtly greater than Bd — >■ D'^D°K~, where the ratio in Eq. ((55|) is the other 
way around. This is due to the contribution from the off-shell part. The off-shell contribution in high momentum 
region favors Dsj{2700) — >■ D*K over Dsj(27Q0) — >■ DK, as one can see from the strong interaction matrix elements 
in Eq. (|3ip. The former coupling is quadratic in momentum, while the latter is only linear. The numerical results 
show that the off-shell effect is about 10%. This correction also echos our assertion that the contribution of Dgj {2700) 
should be treated in a three-body picture. 

The effect of off-shell Ds poles can be read from Scenario I shown in parenthesis. For the first two categories, only 
D* pole contributes, while for the latter two categories, containing the current generated D*K, the Dg pole starts 
to contribute as well. This explains why modes in Category 3 and 4 have larger branching fractions in Scenario I'. 
It is interesting to note that all branching fractions in Scenario I' are deficient in explain experimental results. The 
Z^sj (2700) resonance provides an important source for the non-negligible three-body branching fractions of current- 
produced modes. Comparing with Scenario I, one finds the interference between Z?sj(2700) and Ds* poles are 
not negligible. For example, in the 5° — > D*^D*'^K^ decay rate (see Category 4 in Table W} . the Dl*' and Dgj 
contributions are ~ 0.15% and ^ 0.52%, respectively, while the total predicted rate is ~ 0.94%, which implies a 
fairly effective constructive interference between these poles. If the Dsj width were narrow, we would expect the 
interference effect to be negligible and it would be enough to consider a real Dsj(2700) in two-body final states. 

After the above discussion, one can now understand the total branching fractions in Scenario I by combining 
contributions of three different poles (sec the left portion of Fig. [3]). The contribution of Dgj {2700) dominates over 
Di*'. To first order. Category 2 {D*DK) and 4 {D*D*K) have the same branching fractions from Dsj{2700) and 
are larger than Category 1 {DDK) and 3 {DD* K). Ds' poles further split the two categories that have almost the 
same Dgj {2700) contribution. Consequently, modes in Category 4 {D*D*K) have larger total branching fractions 
than Category 2 {D*DK), and similarly for Category 3 {DD*K) and 1 {DDK). The three different poles form the 
hierarchy of total branching fractions of the four categories in Scenario I. 

Now we consider the situation that both the measurements of BaBar and the contribution of Dsj{2700) measured 
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Scenario I (I ): 
Pole Contribution Only 



Modes with K 



Modes with K* 



Modo(f) Bj{B,^f){%) Br{B,-if)(%) ATf/T,{%) 



Modo(f) Bj{B,^f){%) BT{B,^f){%) ATf/T,{%) 






0.19 ±0.12 ±0.04 0.04 ±0.02 ±0.01 0.17 ± 0.10 ± 0.03 

(0.06±0.03±0.01) (0.03±0.02±0.01) (0.09±0.04±0.02) 

0.19 ±0.12 ±0.04 0.04 ±0.02 ±0.01 0.16 ± 0.09 ± 0.03 

(0.05±0.03±0.01) (0.03±0.02±0.01) (0.08±0.04±0.02) 



D,D°K*' 



(0.07±0.03±0.01) 
(0.06±0.03±0.01) 



(0.03±0.01±0.01) 
(0.03±0.01±0.01) 



(0.08±0.04±0.02) 
(0.08±0.04±0.02) 



DID°K- 
D'D^K 



0.64 ±0.43 ±0.13 0.09 ±0.05 ±0.02 0.38 ± 0.23 ± 0.08 

(0.07±0.03±0.01) (0.06±0.03±0.01) (0.12±0.05±0.03) 

0.62 ±0.42 ±0.13 0.09 ±0.05 ±0.02 0.37 ± 0.22 ± 0.08 

(0.07±0.03±0.01) (0.06±0.03±0.01) (0.11 ±0.05± 0.02) 



DID°K" 
D'D-R'' 



(0.04±0.02±0.01) 
(0.04±0.02±0.01) 



(0.03±0.02±0.01) 
(0.03±0.02±0.01) 



(0.07±0.03±0.01) 
(0.07±0.03±0.02) 



D,D*°K- 
D^D'-R' 



0.30 ±0.18 ±0.06 0.09 ±0.05 ±0.02 0.31 ± 0.21 ± 0.06 

(0.17±0.08±0.04) (0.08 ±0.04 ±0.02) (0.23±0.11 ±0.05) 

0.29 ±0.18 ±0.06 0.09 ±0.04 ±0.02 0.30 ± 0.20 ± 0.06 

(0.17±0.08±0.04) (0.08 ±0.04 ±0.02) (0.22±0.11 ± 0.05) 



D^D*"K'- 



(0.18±0.08±0.04) 
(0.17±0.08±0.04) 



(0.08±0.04±0.02) 
(0.08±0.04±0.02) 



(0.24±0.12±0.05) 
(0.24±0.11±0.05) 



D*D*°K- 
D'D'-R' 



0.89 ±0.59 ±0.18 0.17 ±0.09 ±0.03 0.65 ± 0.39 ± 0.14 

(0.14±0.07±0.03) (0.11±0.05±0.02) (0.23±0.11 ±0.05) 

0.86 ±0.57 ±0.18 0.16 ±0.09 ±0.03 0.64 ± 0.38 ± 0.13 

(0.14±0.06±0.03) (0.10±0.05±0.02) (0.22 ±0.10±0.05) 



DID'"R'^ 



(0.05±0.02±0.01) 
(0.05±0.02±0.01) 



(0.04±0.02±0.01) 
(0.03±0.02±0.01) 



(0.08±0.04±0.02) 
(0.08±0.04±0.02) 



Total 



5.9 ± 3.6 ± l.fl 
(2.6 ± 1.2±0.5£ 



Total 



(1.9±0.9±0.4J3 



TABLE VI. The branching fractions (Sj,r) and width difference (AF/) of the three-body Ds* D^*^ K''^> modes in the scenario 

with only pole contribution. Bj and Br denotes the current-produced decay {Bs — >■ /) and the transitional decay [Bs — >■ /), 

respectively. Z)a,j(2700) is not included in modes with K* in this scenario. The results with only Ds* poles are shown in 

parenthesis. 

^ The contribution from CP conjugate modes is included. 



by Belle are confirnicd in the future. We demonstrate that it is possible to reproduce about all measurements by 
using Scenario II: a pole model with NR contribution in time-like form factors of DK , in addition. Note that the first 
two categories share the same current-produced DK ^ while D* K form factors only appear in Category 3 and 4. Since 
modes in the last two categories already agree with data in Scenario I, using pole model only, no NR contribution is 
introduced in D* K form factors. The branching fractions of modes in the first two categories can be tuned by two 
complex NR parameters in the time-like form factors of DK. These two parameters are fixed by fitting to the observed 
branching fractions of Bd — >■ DdD^K~ and Bd — >■ D'^D'^K" (denoted in the remarks in Table FvT) . The best fit gives 
x^^ = (-75 + 52i)GeV^ and x^f^ = (16 -t- 2i)GeV^ where x'^^ and x^f correspond to the NR contribution in 
DK time-like form factor Fq and Fi , respectively [see Eq. (|38|) ] . Usually the two complex (four real) NR parameters 
cannot be fully determined from two constraints. In this case, however, there is a localized and huge Z?sj (2700) 
resonance contribution in Bd — ^ D^D'-'K^ mode. The NR contribution, which is smooth in phase space, has to cancel 
the I?sj(2700) contribution while maintaining the form factors in other parts of phase space. In other words, the 
phases of the NR parameters are constrained by the complex resonance, while the magnitudes, which control NR 
parts in the off-resonance region, are limited by data. The branching fractions of the fit are shown in Table |Vl where 
100% uncertainties in a; s are included in the first errors. In this scenario, all experimental results, except for the 
explicit disagreement in Bu — ^ DuD^K~ between data, can be explained within uncertainty when NR is included. In 
particular, the Bd -^ D^D^K^ rate is now reduced by a factor of 2 and consistent with data within errors. 



2. Branching Fractions in Bs system and the Width Difference 



After checking the validity of our calculation by comparing to existing data on rates, we move to our main purpose: 
estimating AF^. The relevant diagram is shown in Fig. [3] In Table IVll we show our results in Scenarios I^'-'. Recall 
that bounds on AF^ are related to rates [see Eq. (J15dp ]. The branching fractions of current-produced modes and 
transition modes are also shown, and can be read from Bj{Bs — > /) and B-j-{Bs — ^ /), respectively. For simplicity, 
only modes with K^*'> are shown and the results of modes with iiT^*' can be derived from their CP conjugates. As 
noted before, since CP is conserved in this work, B{Bs -> /) = B{Bs -^ f) and AF/ = ^r^. The total AF//Fs 
contains modes in the table and their CP conjugates, so it is twice the sum of the listed AF^/Fs in the table. 

Before discussing AFg, we first look at branching fractions of these modes. Current produced modes in Bg decays 
are SU(3) related to modes considered previously. Their rates are similar. For example, Bg —>■ D*D*K modes have 



17 



Scenario II: 
Pole contribution + NR in DK tinic-lilce form factors 



Modes with K 



Modes with K' 



Mode(f) B^(B,^f){%) BT{B,^f)(%) AT f /T ,{%) 



Mode(f) Bj{B,^f){%) B-riB, -^ f){%) AT f /T ,(%) 






0.093 
O.OQ'^ 



± 0.02 
±0.02 



0.04 ± 
0.04 ± 



0.02 ±0.01 
0.02 ±0.01 



0.08 ±0.15 ±0.01 
0.07 ±0.13 ±0.01 



D,D°K*^ (0.07±0.03±0.01) 
DsD-R*" (0.06±0.03±0.01) 



(0.03 ± 
(0.03 ± 



0.01±0.01) 
0.01±0.01) 



(0.08±0.04±0.02) 
(0.08±0.04±0.02) 



D^D"K~ 



o.ai: 

0.29 



+0.- 
-0-: 



±0.13 
±0.13 



0.09 ± 
0.09 ± 



0.05 ±0.02 
0.05 ±0.02 



0.11 ±0.38 ±0.02 
0.11 ±0.38 ±0.02 



DID"K'^ (0.04±0.02±0.01) 



(0.03 ± 
(0.03 ± 



0.02±0.01) 
0.02 ±0.01) 



(0.07±0.03±0.01) 
(0.07±0.03±0.02) 






0.30 ± 
0.29 ± 



0.18 ±0.06 
0.18 ±0.06 



0.09 ± 
0.09 ± 



0.05 ±0.02 
0.04 ±0.02 



0.31 ±0.21 ±0.06 
0.30 ±0.20 ±0.06 



DsD'^K'^ (0.18±0.08±0.04) (0.08±0.04±0.02) (0.24±0.12 ±0.05) 
D^D*-K'° (0.17±0.08±0.04) (0.08±0.04±0.02) (0.24±0.11 ± 0.05) 



D'D*°K- 



0.89 ± 
0.86 ± 



0.59 ±0.18 
0.57 ±0.18 



0.17 ± 
0.16 ± 



0.09 ±0.03 
0.09 ±0.03 



0.65 ±0.39 ±0.14 
0.64 ±0.38 ±0.13 



DID'^K'- (0.05±0.02±0.01) (0.04±0.02±0.01) (0.08±0.04±0.02) 
D*Z)*"/?*" (0.05±0.02±0.01) (0.03±0.02±0.01) (0.08±0.04±0.02) 



Total 



4.5±4.4±0.fl 



Total 



(1.9±0.9±0.4J3 



TABLE VII. The branching fractions (Bj^f) and width difference (AF/) of the three-body Dg D'^*^ K^*^ modes in Scenario II 
where DK time-like form factors have NR contribution. The notation is the same as in Table IVTl 
^ The contribution from CP conjugate modes is included. 



largest rates (^ 0.88%) as the Bu^d — > -D* d^*^ modes. However, the transition modes are new. Their rates are 
sub-percent or smaller. Note that while current-produced modes with K are dominated by Dsj(2700), transition 
modes do not change significantly when Z)sj(2700) is included. For instance, without Dgj the branching fraction 
of current-produced mode Bg -^ D*D^K^ drops from 0.64% to 0.07%. In contrast, it drops only from 0.09% to 
0.06% for the branching fraction of transition mode Bg — >■ D*D^K~ . The distinct behavior is not surprising because 
Bs — >■ D|i5s,7(2700) rate (before D^j — )■ D^K~) is relatively suppressed compared with the Bg — >■ D*D* ones (before 
£)* — >■ D^K~) (see Sec. III. A). As we will see later, the different roles played by these poles will be useful to enhance 
AFs through interferences. 

As the branching fractions of transition modes are not tiny, one would expect a non- negligible AF^. The AF//Fs 
of three-body modes range from 0.07% to 0.65% as shown in Table VI. The last two modes with K have the largest 
AF/ as their rates are largest. In this scenario, the total AF^/F^ is 



AF^/Fs(I?i*'**^L>(*'**)) = (10.4 ±2.5 ±2.2)% 

AF,/F^(D(*)i)(*)i? + Di*'^D'-*^K) = (5.9 ± 3.6 ± 1.2)%, 

ATjTs{Di*W^*'^K* + D^*^D^*^K*) = (1.9 ± 0.9 ± 0.4)%, 

AFjFs = (18.2 ±7.0 ±3 



(39) 



Clearly, the AF^ of three-body modes is comparable to two-body modes. The AF^ of three-body modes is mainly 

(^\ — t^\ 
comprised of modes with K. It shows that the approximation in which Dl Dl modes saturate AF^ is dubious. In 

addition, Eq. p9p agrees with the short-distance calculation in Eq. ([5]) within uncertainties. There is no evidence of 
the violation of short-distance result and the underlying OPE assumption. 

The interference between _Dsj(2700) and Ds can be studied by comparing Scenario I with Scenario I' and the 
result of i)sj(2700). The full treatment of modes with K in Scenario I, where _Dsj (2700) and Ds are taken into 
consideration simultaneously, gives AFs/Fg ~ 5.9%. On the other hand, one can treat Z)sj(2700) and D)* separately 
and sum their AFs/F^. The contribution of Ds only (Scenario I') can be read from the Table. For Z)sj(2700), its 
contribution can be estimated from the two-body calculation (see Sec. Ill A) with narrow width approximation. We 
further check that it decreases from the two-body result of 1.9% to 1.7%, when full three-body calculation is imposed. 
In the case that l)sj(2700) and Dg are sum separately, the total AF^/F^ of modes with K is only 2.6% + 1.7% = 4.3%, 
smaller than 5.9% in Scenario I. The difference, which is about the size of the Dsj contribution alone, shows that 
there is considerable interference between _Dsj(2700) and Dg poles. Such interference can be understood as foUowes. 
As depicted in Fig. [31 the D'^*'>K^*^ pairs emitted by the current-produced l)sj(2700) pole interfere with the same 

states from the transited Dg poles in transition diagram. Unlike the highly suppressed Bg — >■ Dsj transitions (see 
Table ITV)) . the Bs — >■ Dg transitions are sizable (see Table HI]) . leading to enhanced Bg — Bg mixing and AF^. In short, 
AFs receives the interference contribution from current-produced l)sj(2700) pole (from Bs decays) and transited Dg 
poles (from Bs decays), which bypass the mismatch of current-produced and transited Dsj in two-body modes. In 
total, diagrams containing Dsj pole contribute more than those with Ds poles only. 



18 



Scenario III: 
-Dsj{2700) is included in all modes 



Modes with K 



Modes with K* 



Mode(f) 



(exp.) 



/)(%) Bt{B, 



/)(%) 



Ar^./r,(%) 

(limit) 



Mode(f) Bj{B, -* /)(%) Bt{B, -^ /)(%) 



Ar^7r,(%) 

(limit) 



D,D°K- 
D,D~K° 



0.09 
0.09 



+0.: 
+0.; 



±0.02 
±0.02 



0.04 ± 
0.04 ± 



0.02 ±0.01 
0.02 ±0.01 



0.08 ±0.15 
0.07 ±0.13 



±0.01 
±0.01 



D,D°K*- 0.10 ± 0.05 ± 0.02 
D.D^K'" 0.09 ±0.05 ±0.02 



0.03 ±0, 
0.03 ±0, 



02 ±0.01 
01 ±0.01 



0.11 ±0.05 
0.10 ±0.05 



±0.02 
±0.02 



D.13 ^ "■-^"' 

ill ±0-^3. 



DID"K- 



0.3i: 
0.29J 



0.09 ± 
0.09 ± 



0.05 ±0.02 
0.05 ±0.02 



0.11 ±0.38 
0.11 ±0.38 



±0.02 
±0.02 



DID"K' 



0.27 ± 
0.25 ± 



0.13 ±0.06 
0.12 ±0.06 



0.06 ±0. 
0.05 ±0, 



03 ±0.01 
02 ±0.01 



0.23 ±0.11 
0.21 ±0.10 



±0.05 
±0.04 



D,D''°K 



DsD'-K" 



0.30 ± 
0.29 ± 



0.18 ±0.06 
0.18 ±0.06 



0.09 ± 
0.09 ± 



0.05 ±0.02 
0.04 ±0.02 



0.31 ±0.21 
0.30 ±0.20 



±0.06 
±0.06 



DsD*"/^*" 0.28 ±0.13 ±0.06 0.10 ± 0.05 ± 0.02 0.32 ± 0.16 ± 0.07 
D^D*"/?*" 0.27 ±0.13 ±0.06 0.09 ± 0.04 ± 0.02 0.31 ± 0.15 ± 0.07 



D'D*°K 



D'D'-K^ 



0.89 ± 
0.86 ± 



0.59 ±0.18 
0.57 ±0.18 



0.17 ± 
0.16 ± 



0.09 ±0.03 
0.09 ±0.03 



0.65 ±0.39 
0.64 ±0.38 



±0.14 
±0.13 



D;D*"ii"" 0.23 ±0.11 ±0.05 0.05 ± 0.03 ± 0.01 0.21 ± 0.10 ± 0.04 
_D*D*"i?*" 0.21 ±0.10 ± 0.04 0.05 ± 0.03 ± 0.01 0.20 ± 0.09 ± 0.04 



Total 



4.5±3.0±0.fl 



Total 



3.4 ± 1.6 ±0. 13 



TABLE VIII. The branching fractions (Bj^r) and width difference (AF/) of the three-body Ds'D^*^K^*^ modes in Scenario 
III, where Dsj (2700) is included in all modes. The notation is the same as in Table IVll 
^ The contribution from CP conjugate modes is included. 

One can bound the width difference in Table IVll bv Eq. ()15dp . For example, the ATf/Ts is bounded to be 0.77% and 
0.08% for D*D*°K~ and D*D*°K*~ modes, respectively. Comparing to AT/, we see that the bounds in modes with 
K are higher within 20%, while they constrain ATf very well for modes with K*. The accuracy of AVf estimation 

in modes with K* has to do with the virtual Ds poles. The pole contribution of Ds is almost real and so are 

the resulting amplitudes. As a result, the suppression from the inequality of Eq. (J15bp is tiny for modes with K*. 

— /^"^ 
This demonstrates that the virtual Ds poles are very efficient to mediate the width difference. On the contrary, the 

on-shell Dsj{2700), which plays an important role in modes with K, generates complex amplitudes and result in the 

suppression of AT/ in these modes. 

The results of Scenario II are shown in Table IVnl Only the first four modes with K are different from Scenario I. 

Note that all transition modes and modes with K* are still the same as in Scenario I, since there is no measurement 

at all to call beyond pole model. One can read from the table that the ATf of the first four modes (modes with 

NR) drop by 50% to 70%. The decrease is caused by the reduction of the branching fractions in current-produced 

modes. Morever, the actual ATf moves away from the upper bound in Eq. (J15d|) when the complex NR contribution 

are included. In this scenario, the total AT^/T^ is 

Ar^/r,(i:>(*'**)z)(*'**)) = (10.4 ± 2.5 ± 2.2)%, 

ATs/TsiDl*W^*^K + D'f^D'^*^K) = (4.5 ± 4.4 ± 0.9)%, 
ATs/Ts{D';*^D^*'^K* + Di*)D^*^K*) = (1.9 ± 0.9 ± 0.4)%, (40) 

ATs/Ts = (16.7 ± 7.8 ± 3.5)%. 

Despite the drop of AT f in modes with NR, the total ATg remains similar to Scenario I because these modes are not 
dominant in AF^ . Most features are similar to the previous case. The effect of three-body modes is still non- negligible. 
It is interesting to see that the central value is more consistent to short-distance calculation. The conclusion remains 
the same as in Scenario I. 



IV. DISCUSSION 



We have seen that -Dsj (2700) is important in modes with K. One expects ZJgj (2700) to be non-negligible in modes 
with K* as well. Even though Dsj{27QQ) is not heavy enough to decay to on-shell D^*''K*, its width is wide and its 
mass is close to the invariant mass threshold of D^*^K*. Unfortunately, there is no information about the coupling 
constants of the effective Lagrangian for Psj{2700) -^ p'^*'>K*. Unlike the on-shell D^j (2700) -^ D^*'^K decay, we 
cannot extract the coupling constant of Dsj(2700) — > D'^*^K* directly from data. Thus, for illustration, we set the 
coupling constants in analogy to the coupling constants of D* to Z)(*)A'*^*) vertices 



9d^jDMK' ~ 9DsjD(*)k\ 



9d'DMk 



9D-DMK 



-) W 0.5 X (JD^jDMR- 



(41) 
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Table IVIIII shows the result in this analogy, which we call Scenario III. The results of modes with K remain the same 
as in Scenario II. 

Comparing with the results in previous scenarios, all branching fractions and AF/ increase. As before, the effect of 
l)sj(2700) is stronger in current-produced modes than in transition modes. In particular, current-produced modes in 
Category 2 {D*DK*) and 4 {D*D*K*) are very sensitive to the appearance of Dsj{2700). Their branching fractions 
rise at least four times. This large effect of current-produced Dsj{2700) in Category 2 and 4 is similar to modes 
with K. If there is a measurement of modes in these two categories, it is possible to extract 9Ds,d('-^K' by fitting 
to branching fractions. The go^jDMK'- ™- return could help the identification of Z)sj(2700). The current-produced 
modes with K* have branching fractions in the order of 10^'^, similar to modes with K. 

The rise of branching fractions in current-produced modes lead to the increase of AF^.. Following the trend of 
branching fractions, AF/ in Category 2 and 4 have significant increase compared with the other two. In this scenario, 
the total AFs/Fs is 

AFs/Fspi*'**)!)^*'**)) = (10.4 ± 2.5 ± 2.2)%, 
l\Ts/Ts{D^*'^D^*'>k + dI*^d'-*^K) = (4.5 ± 4.4 ± 0.9)%, 
AF,/Fs (£)(*) £)(*)/?* + Di*^D^*^K*) = (3.4 ± 1.6 ± 0.7)%, (42) 

ATs/Ts = (18.2 ± 8.5 ± 3.8)%. 

The total AFs induced by modes with K* almost doubles. The effect from three-body modes is strengthen by 
considering the off-shell decay of Dbj{2700) to D^*^K*. For total AF^, the central value returns to the one in 
Scenario I. Total AFs docs not alter significantly as the contribution for modes with K* is not dominant. The result 
still agrees with short-distance calculation. 

The interference in modes with K* is strong. Similar to the discussion in Scenario I, if we leave only D^j (2700) 
and turn off Ds* poles, the resulting AF//Fs of these modes is only 0.3%. It is much smaller than the 1.5% increase 
found in Scenario III (compared to Scenario II). Recalling the result in Scenario I, one finds that modes with K* allow 
more constructive interference than modes with K. For modes with K, the interference is restricted by the on-shell 
Dsj{2700) resonance, which is localized in phase space. On the contrary, the Dgj (2700) resonance becomes off-shell 

and hence smooth in phase space for modes with K*. It is more coherent to the Ds pole contributions and interfere 
with them better. As in the K case, the interference, mediated by the D^*'K* pair, is comparable to the contribution 
of £»sj(2700) itself. 

We show that the branching fractions of these modes are in the order of 10^"^ to 10^^. Recall that there is no 
corresponding measurement in current-produced modes with K* and in all transition modes. For current-produced 
modes with K*, they can be studied in Bu^d system in analogy to modes with K. These branching fractions should 
be measurable with current data collected by the B factories. On the other hand, Bu^d systems have more different 
behaviors in transition modes. B^^d transit to D^*'>'k pairs instead of D'^*'> K'^*\ The D^*'>Tr pairs can be produced 
either from nearly on-shell D* or from other on-shell intermediate resonances. One expects the transition modes in 
Bu^d are enhanced than in Bg. In fact, semileptonic modes with B^^d —^ D^ nu,d transition have been measured '34*. 
The branching fractions are around 0.5%, much larger than the transition modes in this work. For the purpose of 
estimating the width difference, AFg can be bounded by Eq. (J15d[) when current-produced and transition modes are 
measured. Independent of AFs, experimental studies of these modes will be interesting enough in their own right. 

So far we fit the decay constant of Dgj (2700) by its contribution to B^ —> DuD^K~ as measured by Belle. If 
future experiments favor the result of BaBar and lower the contribution of Dsj (2700), then the decay constant will 
be smaller. In such case, the branching fractions of modes in Category 1 and 2 in pole model become smaller and 
may be consistent with experiments without resorting to NR contribution. Nevertheless, the branching fractions of 
modes in Category 3 and 4 will be deficient. Similar to Scenario 2, one can then add NR contribution in the time- like 
form factors of D*K to fit the observed branching fractions. Although there are more NR parameters in D*K form 
factors, one can extract information in the Dalitz plots, especially the interference between the continuum and the 
Dsj{2700) resonance. These can be studied after future measurements are done. 

In principle, modes with ss, such as 77*^'-' , w, and (p, can also contribute to AFg. These modes are difficult to calculate 
because they mix current-produced, transition, and color-suppressed diagrams together. Nonetheless, we find that the 
contribution of these modes are small. The phase space is suppressed and the number of modes are less. We estimate 
the contribution to AFs by Ds'Dsil^'^ modes with color-allowed diagrams only. The effect is less than 0.7%, which 
is negligible. 

We have shown that the effect of three-body modes could be sizable. It is interesting to see if other high-order 
modes could have similar effect on AFs. Note that the phase space is gradually saturated from DgDK mode to 
D*D*K* mode, and the effect of high-order modes may be limited. Fig. U shows the diagrams of possible four-body 
modes. The first type of diagram (left diagram in in Fig. S]) can produce d'^*^ D^*^ K''*^ K'^*\ but the two K mesons 
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B, 



K„ 



Da 



D„: 



K„, 




(a)4-body diagram type I 



(b)4-body diagram type II 



FIG. 4. Left: The first type 4-body diagram. Both of current and transition part have an qq pair. Right: The second type 
4-body diagrams. AU qq pairs hes in current, or in transition for Bs decays to this mode. 

cannot be simultaneously in excited states because of insufficient phase space. The amplitude of this diagram can 
be calculated with the same form factors as in three-body modes. We roughly estimate the branching fraction of 
this type of diagrams, which is two orders of magnitude smaller than three-body modes. Given that the number of 
]j(*) ij{*) Xi*) Xi*) niodes is 48, only 0.5 times more than 3-body diagrams, the contribution of these diagrams are still 
negligible. The second type may involve pions and could have a larger phase space. We calculate the dimensionless 
fraction of phase space area 



1 ^^ (4-body) 
^* (3-body) 



< 10" 



TO, 



(43) 



where .4* is the phase space area. This ratio strongly suggests that the effect of 4-body modes is negligible. Even if 
the branching fractions of current amplitudes are large, the branching fraction of transition diagrams may not be as 
large as in current amplitudes. It should be safe to estimate AF^ up to three-body modes. 



CONCLUSION 



In conclusion, we have estimated the long-distance contribution to AFs of the Bg 



contributions by two-body Di*^!)^*^ modes. The AF^/F^ by these modes is (10.2 ±2.2 ±2.1 



Bg system. First, we revisit the 
^, which decreases from 
previous result in Ref. [21| . More precise measurements in Bs system can help extract more accurate parameters and 
improve the theoretical prediction. After including 15*0(2317), Dsi{2A&Q)^ and Dsi(2536) resonances, the AF^/Fs 
change only slightly to (10.4 ± 2.5 ± 2.2)%. 

For the three-body DsD''*'^K'^*'> modes, factorization formalism with form factors modeled by Dg and Z)sj(2700) 
poles and non-resonant (NR) contributions, if necessary, are used. The branching fractions predicted by pole models 
are consistent with experiment in two of the four categories, while the agreement in the remaining modes with data 

can be achieved by including NR contribution. Three-body modes can bypass some difficulties in two body modes. 

- - 1^\ 

In particular, sizable constructive interference between Dsj and Ds poles, which is impossible for two body modes, 

arc found. 

Our results on AF^. in three scenarios are summarized in Eq. (p9| . (|40)) . and (|42|) . Although the three scenarios 
have different theoretical assumptions, it is of interest to note that the resulting AF^ values are similar. Thus, 
we give the following concluding remarks. First, the total AFs agrees with short-distance calculation. In other 
words, long distance contributions from h — >■ ccs decays do not enhance AF^ (or the real part of Fi2,s) significantly. 
This demonstrates that the short-distance result and the assumption of OPE are reliable. If the anomalous dimuon 
asymmetry with sizable AFs is confirmed in the future, the enhancement in AFs niust have origins from new physics. 

Second, we find that the effect of three-body modes (~ 8%) is comparable to two-body modes (^ 10%). The 
assumption that two-body decays saturate AFs, receives a considerable correction. This correction comes from both 
L>sj(2700) and off-shell D^*^ poles. 

We end our conclusion by pointing out some experimental issues where progress can be made in the near future. 
Two body modes in Bs decays need to be measured with better precisions (see Sec. III. A). For three body modes. 
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F(0) a b 

O.f 

pB,D, 0.67 ±0.03 1.24 0.46 

yB.o; 0.77 ±0.04 1.42 0.68 

Ao'^'' 0.65 ±0.03 1.37 0.63 

Af'^'- 0.62 ±0.03 0.77 0.11 

Af=-°^* 0.59 ±0.03 1.27 0.56 

TABLE IX. The transition form factors for Bs — > D^* used in this work. 

up to now, there is no measurement of transition modes, nor on modes with K* in B^.d system. Even the available 
measurements in current-produced modes with K contain inconsistencies. In particular, the 2.2a difference between 
Belle and BaBar in B^ — > D'^D'^K~ mode has to be resolved. From Tables in Sec. Ill and IV, we see that many modes 
remain to be found or confirmed experimentally. For example, Bs — >■ D*D''*'^K'^*'> rates are predicted at the percent 
level and can be observed soon. Note that the modes with D^*'>K* will be useful to extract the D^j strong coupling. 
Although the measurements of two and three-body decay rates are useful for refining the theoretical prediction and 
to set bound on AFg, these modes are of interest in their own right. We hope that (Super-) B factories and LHCb 
can complete the measurements of these missing modes. 

Note Added. After the completion of this paper, we noticed the work of Ref. 42], which pointed out that penguin 
contributions to B^ — >■ J/ipfp could reduce somewhat the need for enhanced AFg. It also reiterates the point made 
in the second reference of Ref. Q that there is no indication of large or ill-behaved corrections to the short distance 
expansion (or Heavy Quark Expansion). 
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Appendix A: Basic decay constants and form factors 

The value of basic parameters are summarized in this section. We take Wilson coefficients ci = 1.081 and C2 — 
—0.190 with naive factorization. This corresponds to 

ai = 1.02 ±0.10, (Al) 

where we estimate a 10% uncertainty. The decay constants of D^.d and form factors of B^.d ~^ Du,d are given in 
Ref. 0. 

For calculating Bg ^ Dl transition form factors, we use the same method in Ref. [29]. The Dl decay constants 
are taken to be 

/z5, =260±13MeV, 

/i5; ==260±13MeV. ^ ' 

The decay constant of Dg is consistent with the measured values in Ref. 34] . The decay constant of D* should be the 
same as Dg in heavy-quark limit. Using these two decay constants as constraints, we calculate the transition form 
factor, which is parametrized as 



1 - ag2 + 6g4 



^"-"-(^^)= . :i\u.. - (A3) 



The three parameters f (0), a, and b of different form factors are given in Table HXl 
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Appendix B: Some Conversion and Transformation of Form Factors 

Table 1X1 provides the conversion of our notations to the usual notations of standard form factors. 





Ds 


Dso 


Dl 


D^i (2460, 2536) 


f-D(*)s 


Ids 


fosO 


fo's 


-fosi 


pB^-Ds 

pB.-Ds 


pB^Ds 
pB.Ds 


jpB^DsO 
^0 







Ft 


j-D's 


Ft 


,Ti*s 


pB, 


,T>*s 


-'2 




pB. 


,Ti*s 


-'o 





yBsD*s "^B5+"^zj5i aBsD* sO 

aBsD s m£jg— TTi^gi ^/BsD sO 

aBsD*s m^ g-\-m]jsi rrBs D* sQ 

aBsD's r/BsD'sO 

^0 ^0 



TABLE X. The conversion of the form factors notation in this worlc to the usual notation in the literature. 

If CP is conserved, the form factors of current produced particle pair and antiparticle pair can be related. For the 
standard form factors, the transformation reads 

JDs.D's.Dsl = +JDs,D*s,Dsl, JDsQi^DsV = ^JDsO,DsV, 
pBs(D's,Dsl)^^2^ = +Ft^'''''''''\q'), Ft^'^'ig') = -Ft""'"' (q'), 

where Dsi and Dgv are the CP-even and CF-odd states of the linear combination of £'si(2460) and Z?si(2536). The 
relations for form factors in Eq. ((2T|) to Eq. ((23)) are 

Fo7i<l') = -^oW), 

AVP,vy^q^) = -A^-^iq'). 
The transformations for transition form factors from Eq. (j24p to Eq. ((26)) are 

yBsPP^BsVP.BsVV _ _yBs'PP,B,'PP,B,'PP 

— — — (B3) 

aB^PP,B^VP,BsVV ^ ,aB,PP,B,VP,B^VV ^ ' 

Compared with Eq. (|B2p . there is one additional minus sign coming from the pseudoscalar B^ meson. 

Appendix C: Pole Contribution to Form Factors 

For simplicity, we only list the contributions from Dg and I?* poles. The contributions of Dsj(21QQ) have the same 
forms as D*, but with different mass, width, and strong coupling constants. 

In the time-like DK transition form factors, D* is the only possible pole. But there is an ambiguity in the matrix 
element (Dif |i>Ceff l-Dint) when D* goes to offshell. The matrix element is given by 

(i:)(pi3)i4:(px)|i/:cff|Ant(PD*>eD*)> =<^\n\.-{^{pK-VD)+OLq), (CI) 

where a is undetermined since the associated term is zero when Z?* is on-shcU. According to this matrix element, the 
pole contribution to time-like form factor becomes 

vDK,2s gP-DPfPt^^mint* 1 

^1 '-^ -' "^ 2 r~~. r o' 



pPK,„2-^_ gP'Ppfp:,,,rnint* q^- mf„t, , 9^ _ 1 



^0 (1 ) --3 — z:2 — r-— — T^ — y — z::2 \zi2 — z:t'^^o))^ 
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where mint* and Tint* are the mass and width of the D* pole, respectively. If a is nonzero, A^^ {q^) will increase as 
q^ increase. Such energy dependence is unnatural for form factors. We hence set a as zero. Once a is fixed, we have 
the following pole contribution to transition form factors 

V gO'DP N^ ZV '^ 



g2 _ ^2 + imi„t*ri„t* 2 niB,, + mint* 



B, y '"int* 

BsDK „ ,2 

= (3 — —2 r—. T^ ):T73^— 2 3Tl(9 (pd-Pk) -2 OT ) X 



rriB, q^ - rn^nt* + *mint*rint* 2(g'2 + g^ _ rn?^ 



,1^3^ + mint* ,BsD' , /I '^S^ ^ "^int* N ^2" 2mint* . Be-D* ^ 



( '"^- ";;"'"** <°^' + (1 - -^^ pnu ^_^ — taii^A^^^'), (C3) 

g'2 q'^ TUB, + mint* 9^ 

^2 / 9D'DP \~^I I \ aB,D' 

= y— 2 r- T^ )^r{mB, + mint*)Ai 

rriB, q - mf-au. + *mint*l int* 2 

aB^DK 2 ^2 ™2 
= {— 2 -]-■ T^ '^^ 2-[ o 2 (™ge + mint,)^i 

m-B, g"^ - mf^^^ + »mi„t*ri„t* mj^ - mj. zmf^,.^ 

+ {q'iPD-pK)- "'^',"'^ qq')^^ 2^11^] + ^ifl^ 

"^int* '^B, + mint* rUB, mB, 

where q = po +Pk is the total momentum of transitioned mesons, and q' = pg — q is the momentum of weak current. 
Other modes receive contribution from both Ds and D* poles. The time-like form factors of D* K are 



'int* 



mD* + vhk q^ - mf^^^ + irriinu^int* 



Af^iq^) = 0, (C4) 



\P*K/2\ 



gp'ppfpi,. 



g^ - m^j + imintrint 

(*) . - 

where m.int(*) and Fintr*) are the mass and width of the Dg poles respectively. The Bg to D*K transition form 

(*) 
factors induced by Dg poles are given by 



TrB^P'K 

^2 f gP'P'P w , ^ aB,P' 

2 = (^ 2 r~- r )i'^B, + mint*)A;^ 

™B, 9 - "^int* + «™int*rint* 

Vl _ , gP'P'P ^ ^2 

*- .,2 _ ^2 , .-^ -n /■; 



I'l-B, ^ ^ "^mt* + *™int*rint* (jnB, + mint*) 
^^ - i-2 2'"":"'" p )(2mint*)<-"- 

m|j g-" - mfjjj, + imint*i int* 



rB.P'K 

-iq'-m'Lu)^^^, (C5) 

"4, 



aBsP'K ^ otaB^D* 

^3 _ . gP'P'P N ZV 



g2 _ ^2^^_^ _|_ imint*rint* (wb, + mint*) ' 



"B, y '"int* 



xB^_D*_ff 

^1 _ / -gP 'PP ^ ^B,_D 

9 V 2 9 

aBsP'K ^2 ^2 



i-2 2T V X^ 

9 - mfj^j + iminti int 



^0 _ ^"^ Ba ~ "^int w -gP'PP ^ ^B,P 

2 2 '^ 2 2 

'B ~ 9 9 "^1 

2 ^2 aBsP'K 



"T-B, ~ 1 1 "^int + ^^^inti int 



_y-mit,A 



The £)ii'* and D* K form factors are parameterized in the same way. The pole part of the DK* time-like form factors 
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are 



I? 
'int 



__ / -'^ ^' - "int ^11^ N 

TOD + TOK* q^ - m\^^^ + iTOint^Fint* 






^2 " (?') = 0, (C6) 



q^ - mf^^ + tmintrint 
And the transition form factors derived from the pole model are written as 

m^B, 1 ^ "^int* + *™int*rint* 

TrB.DK' , f aB.D' 

^1 I -'^jD'DV -. ^2 



(- 



'B, ' q^ - "^Lt* + *™int*rint* (ws, + ^int*) ' 

' g^ - TO?^t* + *TOint*rint* 



^tB.DK' .r 



-(^'-TO^^ntJ^^^, (C7) 



^BsDK* , f n^iB^D* 

g2 - m?^j, + iTOint*rint* (toB, + TOint*) ' 
_ / '2.9DDV X .B,D 

^ v^ 2 — r~ F — '' 1 ' 

1 - mf^t + «™intrint 



™B, 


Af< 


,_DK* 


"^la 


aB, 

^0 


,DK* 



V -.yuuv N J 

2 _„2 ''L2_^2 , ,-^. .r. , ^^0 






rUg — g^ ' 



Finally, the D*K* time-like form factors from Dg poles are 

Vo'''''\q') , -^.fD'Dvfint , 



{niD' + rriK' Y q^ - mf^^ + irriintTint ' 
Vr'^'iq') 



vr^^q') 



-0, 
= 0, 



(tod* +TOk.)2 



D'K*i„2\_( '2g]j*]j*v'm-int*fint* 
\D'K*t 2\ _ I ^4/D*_D'ym.int*/int* , 



<^*(9^) = ( -2fejB.yTO,„,./i^,. ^^ ^^g^ 

g^ - TOfjjj^ + iTOint*! int* 
aD'K* ( 2\ _ I ~2gj|.j).ymint*/int* s, g ^ 4/j|. ji»yTOint<,/int* ^ g 

01 W y — V — 9 9 '■ rn / 9 ^ \ — 9 9 '• T^ / 9 — 

9 - TOf^j^ + iTOintJ int* ^f^t* 9 ^ "^int* + «™int*i int* ^fj^j 

.d-k*^2a_/ 4/D.D-vTOi„t*/int* wl q^ - m\,. ^ m\. q^ 

^02 \q } - 1-9 — —9 — n?z^ — T^ — Ho ^32 ^ 



,jf2 - TO?jj^ + iTOi„t*rint* 2 2TO?„t^ (tod-+TO;<-.)2 
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And the transition form factors are given by the following three equations. The first part is the form factors from 
vector current 



™3 ^ y'. 



q^ - TO?„^^ + imint*rint* mB, + mint* ' 



"-B, y '"int* 






m\^ q^ - m?j^^ + imint*rint, niB, + mint* ' 

Y^^ = ( — ^f^'^:y )(_^!: — ), (C9) 

m%^ q^ - mf„j^ + imint*! int* ruB, + mint* 

m% q^ - m?„t + imi^t^int ^ 



m% q^ - m?^^ + imintFint ° 



The second part which originate from axial currents are 



— = {— 2 —■ r h(^Bs + "^int*)Al 

rriB, q-^ - mf^^^ + imint,rint* 2 

aBsD'K* Af 

^61 ^ f ^JD'D'V . ^ 

rriB, q^ - ra^^t* + *™int*rint* 

1 2 2 aBsD'K" 



2^ g2 I ^^^ 



^ — — ^ 7; — )(-qpK' + 2 — 

int* + *"^int*i int* '^int* 



"^s,^6o = i— — z:2 — TTZ — T^ — ){-qpK* + —z:2 ) ^ (cio) 



9 ^ "^mt* + «™int*rint* mf 



{-{rriB, + TOint*)^f°"° + 2mint*A^°'° - (g'^ - (m| - mf^t*)) ^; }, 

TUB, + mint* 



iBsD'K* 
-{q''^ - {m%- q^)y 



„'2 /™2 „2nn^61 



aBsD'K* o _p 

= ( ^ 2 r~ T^ ) y^B, + mint*)A-^ 

1 - "^int* + *™int*i int* 

, ^Id'O'V w , ^ .BsD* 

= ( -^ 2 r~ r ) ("^S, + mint* ) ^1 

1 - mfj^t* + *mint*rint* 

'ZgD'D'Vmint*Jint* ■. —^2 



m-Ba 


Af^'^* 


mB, 


aB.D'K' 
^21 


"1b, 


^20 



1^ - m?^t* + *™int*rint* mB, + mint* 



"int* 

- 2.g_D * B • V TOint * /int * 
mB, ' g^ - mf^^^ + iTOint*rint 



- . ^S,B* _ ^ 

{-(tob, + TOint,)Af=-° - (g'^ - (to| - mfnt*)) h 2mint*A^'^ 

TUB, + mint* 
aBsD'K* 
-{q -(™B, -q ))— -3 ' 
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and 



aB^D'K* Af aBsD' 



= ( 



q^ - m?„j^ + imint*rint* mB, + m-i^t* 



'Be y ^ '"int* 



^10 ^ / ^JD'-D'-V . 

niB, q^ - TO?„t* + *"lint*rint* 



Af^* 



{-(mB^+mint*Mf°'''-('?"-("^l -™L*)) ^ + 2mi„t*A^=^'} (CU) 



^3 






B^D'K* 

B, ff - "^hit* + ^"lint*rint* q^ - mint* + «™int*rint* 



-I / Ar>sU 

1 /.B^i'*^* jS^B'K* 



+ — (^^^b + ^^ 1 

\BsD'K' 



2<7^ TUB, triBs 

^00° _ r ^2gi3.£).ymint*/int* ^fD-D-V 



rriB, q^ - tn^nt* + «»7iint*rint* g^ - m?jt* + «™int*rint 



-}x 



int* 



-(g'2-(m|^-9^))^ 



BsD'K* 
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